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Abstract. We consider a random walk on the fully-connected lattice with N sites 
and study the time evolution of the number of distinct sites s visited by the walker 
on a subset with n sites. A record value v is obtained for s at a record time t when 
the walker visits a site of the subset for the first time. The record time t is a partial 
covering time when v < n and a total covering time when v = n. The probability 
distributions for the number of records s, the record value v and the record (covering) 
time t, involving r-Stirling numbers, are obtained using generating function techniques. 
The mean values, variances and skewnesses are deduced from the generating functions. 
In the scaling limit the probability distributions for s and v lead to the same Gaussian 
density. The fluctuations of the record time t are also Gaussian at partial covering, 
when n — V = 0(n). They are distributed according to the type-I Gumbel extreme- 
value distribution at total covering, when v = n. A discrete sequence of generalized 
Gumbel distributions, indexed hy n — v, is obtained at almost total covering, when 
n—v = 0(1). These generalized Gumbel distributions are crossing over to the Gaussian 
distribution when n — v increases. 


Keywords: random walk, fully-connected lattice, visited sites, records, Gumbel distri¬ 
bution 

1. Introduction 

In a recent work [1] (referred to as I) exact results for the statistics of the number of 
distinct sites s visited by a random walker up to time t on the fully-connected lattice with 
N sites were presented. Discrete probability distributions, involving Stirling numbers 
of the second kind, were obtained leading to Gaussian distributions in the continuum, 
scaling limit. The present work is a continuation of I in which we study the statistical 
properties of the records associated with these numbers s. The walk is still taking place 
on the fully-connected lattice with N sites but s is now the number of distinct sites 
visited by the walker, belonging to a subset with n sites (see [2] for a study of this 
problem in dimension d = 1 to 3). A record is established each time the walker visits 
a new site among the n sites of the subset (see hgures 1 and 2(a)). To each record we 
associate a record time t and a record value v whereas s gives the number of records. 
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Figure 1. In this example the fully-connected lattice has = 11 sites among which 
n = 5 (squares) belong to the subset. At t = 0 the walker is outside the lattice. A first 
record occurs at the third step with record time t = 3 and record value v = 1. There 
is no record at the next step, t = 4, since the walker visits this site for the second 
time. The two following steps, corresponding to the visit of new sites belonging to the 
subset, bring two new records with, respectively, record time t = 5 and value v = 2, 
record time t = 6 and value u = 3. At t = 7, which is not a record time, the number 
of records is s = 3. The evolution of s as a function of t is shown in figure 2(a). 

We study the discrete probability distributions associated with the stochastic variables 
s, V and t, and their scaling behaviour. 

The study of records [3-7] has been a subject of renewed interest in the physics 
community during the last years [8-17]. In the visited sites problem the records are not 
standard ones since s is a non-decreasing random variable. A record time corresponding 
to the total covering of a graph with N vertices by a random walk is usually called the 
covering (or cover) time of the graph. In the following we use indifferently the terms 
‘record time’ or ‘covering time’. We shall have to distinguish between the partial covering 
time when the record value v is such that Imi^^n^ooV/n = x < 1 (i.e., n — v = 0(n)), 
the total covering time when x = 1 with v = n and the almost total covering time when 
X = 1 with n — V = 0(1). 

Exact results for the mean covering time have been obtained on a lattice with N 
sites in one dimension {ID) with either periodic or reflecting boundary conditions [18]: 

- 1) (periodic), 

t^^\i) = N{N — 1) + {i — 1){N — i) (reflecting). (1.1) 

In the reflecting case i gives the initial position of the walker. The typical size of the 
domain explored by the random walker at time t grows as y/t which explains the long 
time quadratic growth with N. The same problems have been recently solved in the 
case of a persistent random walk [19]. Exact results have also been obtained in ID for 
the mean partial covering time and the mean random covering time [20] which is the 
time needed to visit a given fraction of the N sites chosen at random. 

On the basis of Monte Carlo simulations, the following expressions were 
conjectured [21] for the mean covering time for S> 1 in higher dimensions: 

Iff = A 2 N\n^ N[1 + 0(1/\nN)] (d = 2), 

Iff = AdN In N[1 + 0(1/ In N)] (d > 3). (1.2) 

Logarithmic corrections can be traced to multiple visits of the same sites and this effect 
is stronger in 2D. 
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Actually the leading contribution for d > 3 has been derived earlier using methods 

of probability theory [22]. The amplitude can be expressed in terms of the probability 

of return to the origin. The conjecture (1-2) was also indirectly confirmed analytically 
for d > 2 [23, 24] with the following values for the amplitudes on hypercubic lattices 
with periodic boundary conditions: 

{ l/7r = 0.318... , d = 2, 

1.516..., d = 3, (1.3) 

1.239..., d = 4. 


In mean-field, or d = oo, the amplitude is simply given by A^o = 1 [22-26]. 

The mean covering time for k visits has been studied in lH through Monte Carlo 
simulations [27]. The same method has been used to evaluate partial and random 
covering times in 2D [28]. 

Some exact results are known for the probability distribution of the total covering 
time. An exact closed-form analytical expressions has been derived for a random walk 
on an arbitrary graph [29]. The expression for the complete graph (or fully-connected 
lattice) is in agreement with a conjecture based on small-size exact enumerations [25]. 
There it was noticed that the probability to have a total covering time equal to its 
minimum value is given by the ratio of the number of Hamiltonian walks Zhw to the 
number of random walks Zrw with N steps. Then using a l/d expansion of ^hw for 
the hypercubic lattice [30, 31] it was shown that the first non-vanishing correction to 
Zhw/^rw is of order 1/d^. 

It has been rigorously shown recently [32] that the fluctuations of the total covering 
time for a random walk on the discrete torus in dimension d > 3 are governed by the 
Gumbel distribution [33, 34] in the scaling limit. This result was earlier conjectured 
in [35]. 

Our main results for random walks on the fully-connected lattice can be summarized 
as follows. The number of records s established up to time t is distributed according 
to t 


*S*7V,n('5, t) 


n- jN — n + t'] 


(1.4) 


In this expression N is the lattice size, n the subset size, n- is a falling factorial 
power ( [36], p 47) and ^ an r-Stirling number of the second kind [37]. 

The record value v for a given record time t is distributed according to 


Vv,n(w,i) 


(n- l)^/Ar-n-ft- 1) 
{N-iy-^XN-n + v-lj 


(1.5) 


whereas the record time t for a given record value v < n is distributed according to 


TN,n{v,t) 


n- j N — n + t — 
N‘\N -n + v -ij 


( 1 . 6 ) 


t Note that s is also the number of distinct sites visited on the subset up to time t. The distribution (1.4) 
was obtained in I for n = N 
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In the scaling limit, indicated by ‘ s.L’ (iV —)> oo, n —)■ oo, t ^ oo or n —)■ oo, with 


hxed ratios n/N f, t/N "'= w or v/n = x), s and v have the same mean values SN,n{t) 
and VN,n(t) such that 


s.l. 


s.l 


S7V,n(^) ^ 


(1.7) 


n n 

The probability distributions t) and t) lead to the same centered Gaussian 

density in the reduced variable 


"5 'Sjv,n(^) s.l. V V]\f^nit') s.l. 

^ = a. 


n 


1/2 


n 


1/2 


The mean partial covering time t]\f^n{v) behaves as: 


s.l. 1/1 \ / ^ 1 1 

—= ~ ln(l — X) (x < 1) . 


( 1 . 8 ) 


(1.9) 


The probability distribution T;v,n(n,t) leads to a centered Gaussian density in the 
reduced variable 


]VV2 “ ^ 


(x < 1) , 


( 1 . 10 ) 


at partial covering. The scaling behaviour is different at total covering. The mean value 
of the covering time is given by 

tAr,n(t’) - A^(lnn + 7 ) (u = n > 1), (1.11) 

where 7 = 0.577 215 665... is the Euler-Mascheroni constant. The scaling limit of 
NTN,n{v,t) is the type-I Gumbel extreme value distribution in the reduced variable 


t s.l. / 

- = T 

N 


( 1 . 12 ) 


The crossover from Gumbel to Gauss, in the vicinity of total covering when u = n — v = 
0(1), occurs via a discrete sequence of generalized Gumbel distributions [38], indexed 
by the deviation u from total covering. 

The outline of the paper is as follows. In section 2 the discrete probability 
distributions for the record numbers, record values and record times are obtained using 
generating functions techniques. Their moments are calculated in section 3, leading to 
the mean values, variances and skewnesses. The scaling limit is studied in section 4. The 
results are discussed in section 5. Detailed calculations are presented in five appendices. 


2. Discrete probability distributions associated with records 

We consider a random walk on the fully-connected lattice with N sites (see figure 1 
in I). The random steps are towards anyone of the N sites, with probability 1/N. The 
walker is outside the lattice at t = 0 and the first random step on the lattice takes place 
at t = 1. We study the properties of the number s of distinct sites visited up to time t 
on a fixed subset of n sites arbitrarily chosen among the N. Since s is a non-decreasing 
function of time, a new record for s is established each time a new site is visited on 
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Figure 2. A time evolution of the number of distinct sites visited is sketched in (a). A 
small circle corresponds either to a step outside the subset of n sites or on a previously 
visited site inside the subset. Bigger circles correspond to a record when the walker 
visits a new site inside the subset of n sites. The sequence shown ends with a record 
with value v at time t. The diagrams in (b) and (c) give the contributions to the 
generating function of the probability distribution JN,n{k,t) for the record lifetimes 
in (2.1) and (2.2), respectively. 


the subset of n sites. Thus the value of s at time t also gives the number of records 
established up to time t. A record is characterized by its value v and its time t, i.e., v 
is the record value taken by s at the record time t (see hgures 1 and 2(a)). 


2.1. Generating functions 


In this section we are looking for the generating functions for the probability 
distributions and S]\f^n{s,t). The hrst one gives the probability to establish a 

new record with value v at time t while the second gives the probability to have visited 
s distinct sites (established s records) up to time t. 

Let us hrst write the ordinary generating function for the probability distribution 
JN,n{s, 1) of the lifetime I when the record number is s: 


Jn,i 


s, z) = 


^ ^ k) Z Z 


l=l 


n 


N 


kk’N,n(^^i z') ■ 


( 2 . 1 ) 


In this expression, the hrst two factors give the contribution of the record-breaking hnal 
step (bigger circles) in the diagrams of hgure 2(b), occurring with probability {n — s)/N. 
The last one is the generating function which corresponds to the diagrams of hgure 2(c) 

i 


z') ^ ^ 


1=0 


N — n-\- s 
N 


N 


z = 


N - z{N -n + s)' 


( 2 , 2 ) 


and snms the contributions of seqnences of I steps, each with weight z (smaller circles), 
for which the number of visited sites in the subset remains equal to s, with probability 
[N — n + s)/N. Inserting (2.2) into (2.1), one obtains: 

z{n — s) 


fJ (V,n('S) ■^) 


(2,3) 


N — z{N — n + s)' 

Any history in the (s, t)-plane, ending with a record at (n, t) as shown in hgure 2(a), 
can be decomposed into subsequences of steps at level s ending with a record at s -|- 1 
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as shown in figure 2(b), with s increasing from 0 to n — 1. The record time t is the 
sum of the lifetimes I of the previous records. These are distributed according to the 
different distributions JN,n{s, 1) for s = 0 to n — 1. Thus the record times t are the sum 
of V independent, differently distributed random variables. 

All the possible histories are taken into account if the generating function for 
TN,n{v,t) is written as the product of the generating functions for the lifetimes of the 
successive records, JN,n{.s, 2 ;), for s = 0 to n — 1: 


CXD 

TN,n{v, Z) TN,n{v, H = 


v—1 


n- {z/Ny 


t=i 


n:;J |i - (Viv)(iv - n+s)] 


s=0 


(2.4) 


In the final expression n- = n{n — 1) • ■ • (n — n + 1) is a falling factorial power [36]. 

For SN,n{s, t) the sequence of steps may or may not end on a record. This is taken 
into account in the generating function through a multiplication by )CN,n{v,z) giving: 


00 

<SN,n{^, ^) = SN,n{s, t) = TN,n{s, z) z) 

t=0 


n-{z/Ny 


( 2 . 6 ) 


2.2. Probability distributions 


The probability distribution of the lifetime I of record number s is the coefficient of z^ 
in the expansion of jN,n{,s,z) given by (2.3): 


JN,n(.^y 0 



(2.6) 


This is the geometric distribution with parameter p = 

The expressions obtained for TN,n{v,z) and SN,nis,z) are closely related to the 
generating function of the r-Stirling numbers of the second kind (see [37], equation (25)) 


00 


E 


r + m) ^ 


z 


m 


n 


ra 

k=f) 


[1 


— z{r + k)] ’ 


r > 0, m > 0 . 


(2.7) 


These numbers can be expressed in terms of ordinary Stirling numbers of the second 
kind [36] through ( [37], equation (32)) 



( 2 . 8 ) 


SO that 


f r + / 
\r + m 


r=0 



(2.9) 


Comparing equations (2.4) and (2.5) to (2.7), one obtains the joint probability 
distribution 


TN,n{v,t) 


n- j N — n + t — 1'] 
iV* \n — n + V — l] 


N-n 


( 2 . 10 ) 
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for the record value v and the record time t and the probability distribution 


*S*7V,n('5, t) 


n- jN — n + t'] 

iv\]V-n + 4„_„ 


( 2 . 11 ) 


for the number of distinct sites s visited up to time t on the n-site subset of the fully- 
connected lattice with a total of N sites. When n = iV in (2.11) equation (2.8) of I is 
recovered. 

The last probability distribution satisfies the master equation 

/V — 77-1-9 77 — 9-1-1 

SnAs, t) = - t-l) + -^^W,n(s - 1 , t - 1 ) , ( 2 . 12 ) 

for t > 0, with S'Ar,n(s,0) = and SnA^ < 0,t) = 0. The hrst term on the right, for 
which s remains unchanged, corresponds to steps outside the subset of n sites or on one 
of the s sites already visited in the subset. With the second term the number of visited 
sites increases from s — 1 to s for steps on one of the n — s + 1 sites not yet visited in 
the subset of n sites at time t — 1. 

Let us now look at the properties of Tj^A'^i t)- The probability distribution vanishes 
when V > n due to the falling factorial and when v > t due to the r-Stirling number. 
When summed over t, it gives the probability to establish a record with value v so that: 


CX) 

^TAr,„(n,t) = TNA^yAz=l 

t=l 


1 if n < n 
0 otherwise 


(2.13) 


It follows that TnAv , t) is the properly normalized probability distribution for the 
records times f at a given value of n < n. When summed over v, it gives the probability 
to establish a new record at time t > 0 : 


BnA'^) ~ • (2-14) 

V=\ 


Thus the normalized probability distribution for the value n of a record established at 
time t reads: 


14r,n(w,i) = 




'-N,n 


[V,t) 


(2.15) 


YA=1 TnA'^1 Bn At) 

To close this section let us look for the expression of the joint probability distribution 
in the simple case where v = 1. Then (2.10) leads to 


TnA^^ t) 


n fiV — — 11 

iW\ N-n jj. 


(2.16) 


Making use of (2.8) with {q} = 6kfl gives = r'' and one obtains: 

( 72 \ 72 

^ ~ iv) Jf = J^At^^t) ■ (2-17) 

This is the probability to have the first t — 1 steps outside the subset followed by a visit 
to one of the n sites of the subset at the record time t. 
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3. Moments of the probability distributions 


3.1. Moments of JN,n{s,l) 


Using the expression of the geometric distribution in (2.6) a direct calculation leads to 
the hrst two moments: 

_ °° AT _ °° / N N 

lN.n{s) = TlJ^,As,l) = -. /i,,„(s) = V/0„,„(s,/)=2 --. (3.1) 

n — s ’ \n — s J n — s 

The mean value increases with s and is equal to N when s = n — 1, i.e., for the last 
lifetime before total covering of the subset. The variance is given by: 

^ N 

(3.2) 


It is equal to N{N — 1) when s = n — 1. 


N 


n — s 


n — s 


3.2. Moments of SN,n{s,t) 

Let us start with the number of distinct sites s visited (number of records) up to time 
f > 0. The moments of SN,n{s, t) can be deduced from the bivariate generating function 
obtained in appendix A: 

OO i 


SnAv. ^) = E E [1 + !/(e'/" - 1)]" . (3.3) 

s=0 t=0 

The mean number of distinct sites visited on the subset of n sites up to time t is given 
by the coefficient of /t\ in the ^/-derivative of SN,n{,y, z) ai y = 1: 

/ ,V - , 1 

= n 1 — 

y=l 


SN,n{t) 


dS 


'N,n 


dy 


N-1 

N 


(3.4) 


It is the product of the number of sites in the subset by the probability that a given site 
has been visited up to time t, 1 — (1 — l/Nf. The second moment reads: 


^N.n 


{t) = 


'U' 

d 

dSN,n 1 

.t'-. 

dy 

^ dy \ 


= n^ — n{2n — 1 ) 


y=l 

N-1 


+ n{n — 1 ) 


N-2 


N J ' ' \ N 

Combining these results, the following expression is obtained for the variance: 


= n 


N-1 


-|- n{n — 1 ) 


N-2 


— n 


N-1 


2t 


(3.5) 


(3.6) 


N / ' ' \ N / \ N 

Note that the variance is not proportional to n, which indicates a correlation effect. The 
visits of the different sites of the subset are not independent: when one site is visited in 
a given step, the others are evidently not. 
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The third moment is given by: 


^N,n 


it) = 


A'- 

d 

_t\_ 

dy 


d f dSM,n\\ 3 ro 2 . , 1 fN-l\ 

(--!)+-] [-1^) 


3n{n — 1)" 


N-2 


— n{n — l){n — 2 ) 


N-3 


(3.7) 


N J \ N 

It is needed to evaluate the skewness 71 measuring the asymmetry of the distribution. 
For a random variable X the skewness is the third centered moment, normalized by the 
standard deviation to the third power [40]: 

(X -Xf JO -3XAX^ -X^ , , 

7i = ^ =- , -. 3.8 

(AX2)3/2 (AX2)3/2 ^ ^ 

The skewness of SN,n{,v,t) follows from previous results in equations (3.4), (3.6) 
and (3.7): 

(S) _ 1 -pi + 3np'l-3{n-l)p2-{n-l){n-2)ps + 3n{n-l)pip2-2n‘^pl 


jrlnit) = 


ni/2 


[pi + (n - l)p2 - npl 


where pt = 


N -k 
N 


(3,9) 


3.3. Moments of VN^n{v,t) 

In order to study the moments of t) let us hrst look closer at the relations between 

7jv,n, VN,n and S^^n- First, the probability to establish a new record with record value v 
and record time t is also the probability to have visited u — 1 distinct sites at time t — 1 
and to visit a new site among the n — v + 1 remaining ones at time f, thus we have 

TnA'^^ = SN,n{v - 1, t - 1) - -= *S'jv,n(n, t) “ ^ SN,n{v, t - 1) , (3.10) 

where the master equation (2.12) has been used. Inserting the last expression in (2.14) 
leads to 

BN,n{t) = ^ SN,n{v, t) -^ ^V,n(n, ^ “ 7) “ jy ^ >^V,n(n, t - 1) 

V=1 V=1 V=1 

n AT _ ^ 1 

= t) - ^ ^ - 1) ^ V ^7V,n(n, t - 1) 


?;=0 


i ;=0 


v=l 


1 


1 


n SN,n{t — l) n —1 
A ~ ' 


t-i 


(3,11) 


N N \ N 

when f > 0. The contribution of the terms u = 0 added in the two hrst sums of the 
second line vanishes since SN,n{^-it) = — 1). The hnal expression, which 

follows from (3.4), is the product of the probability to visit a site belonging to the 
subset at t, by the probability that it has never been visited before. Making use of this 
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expression in (2.15) finally gives: 

T. / N N f N 

VN,n[vR) = 


t-1 


n-jN — n + t — 1 


n \N — 1J \N — n + V — 1 j 

{n — 1)^^ r N — n + t — 1 


iv-ilvr ^ II t-1). (3.12) 

{N-ly ^ [N-n + v - 1) 

Since S'Ar,n(s,0) = (5s,o, Ihe last equation gives VN,n{v, 1) = (5„,i, as expected. 

According to (3.12), V^^nivyi) is given by S' 7 v-i,n-i(s,t — 1) at s = u — 1. Thus 
their mean values are shifted: 


Wv,n(t) = SN-l,n-l{t - 1) + 1 . 


(3.13) 


The variance and the skewness, which are both functions of n — VN,n{t) are not affected 
by the shift. They are simply given by: 

^^V,n(t) = - 1) , 7i^^,n(t) = 7S-l,n-l(t “1) (t > 1) • (3.14) 

Note that, since VN,n{v, 1) = the variance vanishes and the skewness remains 
undehned when t = 1. 

3 . 4 . Moments of T]\i^n{v,t) 

The hrst derivative at z = 1 of the generating function in (2.4) gives the mean value of 
the record times (covering times) as a function of the record value v: 


^JV,n('y) = y^jTN,n{v,t) 


d% 


N,n 


t=l 


% 


v—1 




dz 

N 


Tn.u dlnTiy 


N,n 


z=l 


2 = 1 


7V,n 


^ ^ ~ — n + k) 


z dlnz 
= N{H^ - H^_,) . 


(3.16) 


2 = 1 


Here Hn = Y12=i 1/^ is a harmonic number with Ho = 0 [36,41]. The mean record time 
is the sum ^v,n(s) of the mean lifetimes in (3.1). The second moment is given by 

V—1 

X] M - 


72-TT + 2 rp f d f OTn^tI 

TN,n{v,t) = —{Z 


t=l 

% 


N,n 


dz 

N 


dz 


2=1 


1 . 

dz 


N 


^ N — z{N — n + k) 


2=1 




'ii-i 

N — z{N — n + k) 


V — 1 


^ ^ N{N-n + k) 

+ /iV,n 2_^ \AT _ -yf AT _ ^ mi2 


2=1 


i^jN-z{N-n + k)y 


2=1 


= [{Hn - Hn-,f + Hn,2 - Hn-,, 2 ] " N{Hn - Hn-n) , (3.16) 

where Hn,m = ^ generalized harmonic number with i?o,m = 0 [36,41]. 

Thus the variance takes the following form: 

= N\Hn,2 - Hn-,,2) " N{Hn - Hn-,). (3.17) 


This is the sum ni'^) variances of the lifetimes in (3.2), as expected for 

a sum of independent random variables. For the third moment we have: 


4,n(^) = '^t^TM,n{v,t) = 


t=l 


d 

\ d 

1- 

b 

dz 

^dz 

V'^ dz J\ 


(3.18) 


2 = 1 
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Figure 3. The main panels give the time evolution for different lattice sizes of (a) 
the mean value VN,n{t), (b) the variance Avff ^(t) and (c) the skewness ^ (t) of 

the record value v for the number of distinct sites visited by the random walk on the 
n-site sublattice. Full lines correspond to / = n/N = 1 and dashed lines to / = 11^/2. 
Sizes are increasing from bottom to top or left to right in (c). The insets show the 
data collapse obtained with scaled variables (symbols) as well as the behaviour in the 
scaling limit (full line) as explained in section 4. Blue symbols correspond to / = 1 
and yellow symbols to / = The sizes are the same as on the main panels with 

N — n given by 64 — 64 and 64 — 45 (diamond), 128 — 128 and 128 — 91 (square), 
256 - 256 and 256 - 181 (circle). 
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tN,[{Hn - Hn-vf + KHn " h^n-.)(h^n,2 " i/n-.,2) + 

- 3iV2 [{H^ - + {Hn,2 - + N{Hn - Hn-.) . (3.19) 

Inserting these results in the expression of the skewness (3.8), one obtains: 


(T) / N 1 ‘^H‘^{Hn,3 — Hn-v,3) ~ ^H{Hn 2 — Hn-v,2) + Hn — Hn-. 

7l7V,n(^) = 


Ari/2 


[N{Hn,2 - Hn-.,2) - Hn + Hn-.f^ 


(3.20) 


4. Moments and probability densities in the scaling limit 

In the scaling limit —?■ oo, n —)> oo, t —?■ oo or n, s oo, the appropriate combinations 

of variables are $ 


N ~ ^ ' 


t s.l. 

= W, 


N 


V s.l. 

- = X, 

n 


■S s.l. / 

— = X . 

n 


(4.1) 


In this section we look for the behaviour of the probability distributions in this limit. 


4.1. Records lifetime 

Using scaled variables, the mean value of the lifetime in (3.1) and its variance in (3.2) 
take the following form: 

2 1 

(4.2) 


S.l. 


lN,n{s) = 


= = 


J{l-x') 


f{l-x') 


/(I - x') ’ 

Note that they are diverging as total covering is approached, i.e., when x' —?■ 1. The 
probability distribution of the lifetimes in (2.6), which is exponential in /, takes the 
following form 

s.l. 


^7v,n(s,/) ~/(I 
close to total covering. 


(4.3) 


4-2. Number of records and records values 

According to (3.12) the probability distributions *S'Ar^„(s,t) and Vp^^ni^yt) have the same 
scaling limit. There hrst moment follows from (3.4) and (3.13) and scales as n so that: 


SN,n{t) s_J. VN,n{t) sj_. ^ 


W = t/N . 


(4.4) 


n n 

The mean value grows initially as ft (inhnite system behaviour). The approach to the 
saturation value, n, is exponential, with a relaxation time equal to N (see hgure 3(a)). 
The inset shows a good data collapse on the scaling function (4.4). 

The variance given, by (3.6) and (3.14), scales as n too and behaves as: 


s.l. s.l. , / ^ n -.u i -f \ - 2 w 

= nf,w)=e -{l + fw)e 


n 


n 


(4.5) 


f See I for a discussion of the finite-size scaling behaviour on the fully-connected lattice 
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Although the variance is extensive in the scaling limit, we show in appendix B that 
its value per site remains modified by correlations. The initial growth is linear when 
/ < 1, quadratic when / = 1 and the long-time decay is exponential (see figure 3(b)). 
The fiuctuations are at their maximum for t close to N (tCmax = 1.256431... for / = 1 
and WaisLx = 0.991353... for / = l/\/2). The inset shows the data collapse on two 
/-dependent scaling functions given by (4.5). 

For the skewness in (3.9) and (3.14), vanishing as one obtains: 


1/2 s.i. -e-^' +3(1 + fw)e-^^ - (2 +6fw - 2pw + 3pw^)e-^^ 

n ' 7l7V,n(^) = r X N 9 l3/2 

[e”"^ — (1 + fw)e~^^] ' 

When tc 1 the leading contribution to (4.6) is such that 

(4.6) 

and 

1/2 sU 2/1 (f<l) 

7l7V.nRJ _ y)]l/2 ^ 

(4.7) 


«■+"+(*) = ~ (/ = !)■ 

(4.8) 

When w 

/:§> 1 an expansion in powers of e”"' gives: 



l-|(l + /«,)e-"' + 0(e--) . 

(4.9) 


The skewness in figure 3(c) goes through a negative maximum value for t close to n. 
It is diverging in the limits tc —)■ 0 and w ^ oo but this does not mean an increase 
of the asymmetry. On the contrary a closer look at (4.6) shows that the numerator, 
given by the third centered moment and measuring the asymmetry of the distribution, 
is vanishing as w when / < 1 and when / = 1 for tc <C 1 and as e~'^ for tc S> 1. 
Thus the divergences are dne to the cube of the standard deviation in the denominator, 
which vanishes more rapidly in these limits. The inset shows the data collapse on the 
two scaling functions following from (4.6). 

Let us now determine the common scaling limit of the probability distribntions 
SN,n{s,t) and Vvn(n,f). We shall do it starting from the master equation (2.12). The 
scaling behaviour of the variance in (4.5) suggests the introduction of the dimensionless 
variable: 


SN,n{t) s.l. 

= a = 


_ ^1/2(1 _e- 


ni/2 ^ nV2 ^ 

A properly normalized probability density &(a,f,w) is given by n^^^SN,n{s,t) in the 
scaling limit. Then as shown in appendix C, to leading order in an expansion in powers 
of the master equation (2.12) leads to the following partial differential equation 

for the probability density: 

d& e""' d& 

-{l-fe-n^ + cT^ + e- ( 4 . 11 ) 


dw 2 

Rewriting the probability density as 


®(ff,/, K>) = iPlff, ?'(/,!»)] , 


(4.12) 
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d^G _ 

(9cr2 (9 (t2 


where 4>{f,w) is the scaling function defined in (4.5) one has 

dw ’ dw d(j) dw ’ da da ' 
and the partial differential equation (4.11) can be rewritten as: 

dw\d(j) 2 da"^ ) da"^ ^ da 
It is easy to verify that with the Gaussian density 

p-o'^/(20) 

^(^» = , 0 (/,'«^) = - (1 +/w)e 

which is a solution of the diffusion equation on the left, the right-hand side vanishes 
too. 

Starting from the master equation 

JSJ — 77_L7; — 1 79 — 7;_L1 

VnAv, t) = -i - 1) + ,, T VnAv -1,^-1), t>l, (4.16) 


—2w 


,(4.13) 


(4,14) 


(4.15) 


N - 1 


N - 1 


which follows from (2.12) and (3.12), the same Gaussian solution (4.15) is obtained for 
the probability density of the records values with s replaced by v in the expression (4.10) 
of the scaling variable a. 

The scaling behaviour is shown in figure 4(a) for / = 1 and / = l/\/2 at three 
values of ta. A mean-field calculation is given in appendix B. 


4-3. Records times at partial covering 

Let us first consider the case of records times corresponding to a partial covering of the 
sub lattice, x < 1 {n — v = 0{n)). For large values of k the harmonic numbers can be 
written as 

Hk = \nk + j + 0{l/k), (4.17) 

where 7 is the Euler-Mascheroni constant. Thus in the scaling limit we have 

Hn — Hn-v = Inn — ln(n — n) = — In ^1 — — j = — ln(l — x) , (4.18) 

and the behaviour of the partial covering time follows from (3.15) with: 

— ln(l — x) (x < 1). (4.19) 

Note that this expression leads to 

, (4.20) 

s.L 

to be compared to (4.4). The mean value t^A'^)^ shown in figure 5(a), initially grows as 
v/f. The inset illustrates the collapse of the finite-size data on a single scaling function. 


V s.l. , 

— = 1 — exp 
n 


tNA'^) 


N 
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Figure 4. Data collapse obtained for (a) the scaled probability distribution 
n as a function of (n — VN,n{t))/n^^'^ =' a for different values of 

t/N =■ w, (b) the scaled probability distribution as a function of 

(t — tN,n(v))/N^t‘^ =■ T for different values of u/ n =' X and (c) the scaled probability 
distribution NT^^niv, t) as a function of {t—tN,n{v))/N = t' . Blue symbols correspond 
to / = 1 and yellow symbols to / = \l\/2. The sizes N — n are given by 256 — 256 
and 256 — 181 (circle), 512 — 512 and 512 — 362 (triangle up), 1024 — 1024 and 
1024 — 724 (triangle down). The full lines give the probability densities obtained 
in the scaling limit: (a) the Gaussian density ^[a, (j){f,w)] given by (4.15), (b) the 
Gaussian density 0[r, x(/, a;)] given by (4.32) and (c) the type-1 Gumbel distribution 
%iT') given by (4.42). 
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V 



V 



Figure 5. For different lattice sizes, the main panels give the evolution with v of 
(a) the mean value tN,n{v), (b) the variance and (c) the skewness 

of the record time t needed by the random walk to visit v distinct sites on the n- 
site sublattice. Full lines correspond to / = n/N = 1 and dashed lines to / = 

Sizes are increasing from left to right. The dashed-dotted line indicates the value of the 
skewness for n = n in the scaling limit (see equation (4.36) with u=0). The insets show 
the data collapse obtained with scaled variables (symbols) as well as the behaviour in 
the scaling limit (full line) as explained in section 4. Blue symbols correspond to / = 1 
and yellow symbols to / = The sizes are the same as on the main panels with 

N — n given by 64 — 64 and 64 — 45 (diamond), 128 — 128 and 128 — 91 (square), 
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Differences involving generalized harmonic numbers in (3.17) and (3.20) can be 
evaluated in the scaling limit as follows: 


Hr, 


Hn—v,m ^ ^ 


k=n—v-\-l 


1 s.l. 

r du 

1 

1 

1 

_1 

km 

L-v 

1 

1 

1 

1 

_1 

^m-1 


Using (4.18) and (4.21) in the expression of the variance (3.17), one obtains: 

X 


s.l. ., . 

= XU^x) = 


N 


f{^-x) 


+ ln(l — x) (x < 1). 


(4.21) 


(4.22) 


The initial growth of the variance, shown in figure 5(b), is linear when / < 1 and 
quadratic when / = 1. The hnite-size data collapse on the scaling functions is illustrated 
in the inset. 

When X < 1 the skewness in (3.20) behaves as: 


N''h7L(v)= 


3x 


x{2 — x) 

_/2(l-x)2 /(1-x) 


ln(l —x) 


X 




+ ln(l—x) 


3/2 


(4.23) 


When X —)■ 0, to leading order, the last equation leads to 




. 1 . 2xf - 3x/ + X s.l. 2 - / 


(x/-i - x) 


3/2 


[i/(i - /))'/= 


(/ < 1), (4.24) 


whereas 


^ (/ = 1). (4.25) 

Jb 

When X —)■ 1, to leading order, one obtains: 

v,n(^) = ((C < 1) • (4.26) 

The skewness is shown as a function of v in figure 5(c). The inset illustrates the data 
collapse on the two scaling functions following from (4.23). The scaled skewness is 
diverging in both limits but for different reasons. The numerator and the denominator 
of (4.23) vanish when x —?■ 0 and diverge when x ^ 1. Thus the divergence is governed 
by the vanishing variance when x —?■ 0 and by the third centered moment (a measure of 
the asymmetry) when x ^ 1. 

Let us now consider the behaviour of the probability distribution T]si^ri{v,t) in the 
scaling limit. Using (3.10) in (2.12) one obtains the following master equation: 

TnAv, t) = -^-T^,„(u, t - 1) +-- 1, t - 1), t>l, (4.27) 

Taking into account the scaling behaviour of the variance in (4.22) the following 
dimensionless variable is appropriate: 


t — tlSf,n{v) ^ _ t 
771/2 - !■ - jYi/2 


+ ln(l — x) . 


(4.28) 
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Thus the normalized probability density T(r, /, x) is obtained as the continuum limit of 
N^^'^T]si^n{v,t). As shown in appendix D, to leading order in an expansion in powers of 
it satishes the relatively simple partial differential equation: 


^ _ 1 - /(I - x) d‘^‘1 
dx 2/(1 — xY 


(4.29) 


Noticing that x(/, x) in (4.22) is such that 
^ _ 1 - /(I - x) 
dx 


(4.30) 


/(l-x )2 

one can use the change of variables ^{Tyfyx) = n[r, x(/, x)] to further simplify (4.29). 
Doing so, the diffusion equation is hnally obtained: 

d£l 1 


dx 2 dr"^ 

Thus the probability density is Gaussian in the scaling limit and given by: 

.-rV{2x) 


n(r,x) = 




Xif^x) = 


X 


f{^-x) 


+ ln(l — x) . 


(4.31) 


(4.32) 


Note that x(/,x), the scaled variance Atj^^{v)/N in the scaling limit, is diverging at 
X = 1. This means that a new scaling variable is needed at almost total and total 
covering. The hnite-size data collapse on the Gaussian densities is shown in hgure 4(b). 


4-4■ Records times at almost total and total covering 

Let us now study the vicinity of total covering, i.e., the records times at v such that the 
number of unvisited sites u = n — v = 0(1). Then u/n^= 1 — x = 0 in the scaling limit, 
hence the scaling functions in (4.19), (4.22) and (4.23) diverge: the scaling behaviour is 
anomalous. When n is large, inserting (4.17) into (3.15) gives the mean value: 

(4.33) 


u = n — V . 


tN,ni^) - A^(lnn + 7 - Hu) , 
One may write 


n 


l—m 


k=l,n 


m- 


(m > 1), (4.34) 


where C(^) = 1/^™ is the Riemann zeta function. Making use of (4.17) and (4.34) 

with m = 2 in (3.17), one obtains the variance: 


Atl 


N,n\ 


s.l. 


Ar 2 


= ({2)-Hu,2. 


(4.35) 


Thus, when u = n — v = 0{1) the variance is independent of / in the scaling limit and 
scales as instead of N for x < 1. The behaviour of the skewness in (3.20) is obtained 
using (4.17) and (4.34) and reads: 

(D 7i- ^[C(3) ~ -^^, 3 ] 




[c(2) - Hu,2r/^ 


(4.36) 
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Figure 6. Evolution of the covering time probability distribution in the immediate 
vicinity of total covering. In the scaling limit (thick lines) the standard type-I Gumbel 
distribution {u = n — v = 0) is approached via a sequence of generalized Gumbel 
distributions, indexed by the deviation from total covering u = 4, 2,1. Thinner lines 
show finite-size effects for N = 64, 256 and 1024 from top to bottom. 


According to (4.35), the appropriate scaling variable for the probability density is 


now 


-:- ^ ^ — -In n — 7 + hr,, , 

N N n, 


u = n — V . 


(4.37) 


As a consequence, the probability density T(,(r') is given by the scaling limit 
of 7VTy^„(n — u,t) which, as shown in appendix E using the properties of Stirling 
numbers, can be written as §: 


NTN,n{n -u,t) = 


n\ 


n—u—1 


In the scaling limit one may write 

t-i 


u\{n — u — 1)\ -p 
one may 
I -\- u \ 


E (- 1 )' 


n — u — 1 


1 - 


I -\- u -\-1 
N 


t-i 


N 


~ e 






n 




where the last equality follows from (4.37). To leading order in n, one has 


n\ 


n — u — 1 


(n — M — 1)! \ I 
Thus (4.38) leads to 


n\ 


n 




{n-l-u-l)\l\ 


T'(r') = 


0-(“+l)(r'+7--H'u) Q-lR'+l-Hu) 


U\ 


E(-i)' 


1=0 


= -rexp 

u\ 


-{u + l)(r'+ 7 - hr„) - 


(4.38) 


(4.39) 


(4.40) 


(4.41) 


§ Note that when n = N equation (4.38) with u = 0 is in agreement with the probability distribution 
obtained in [29] for the complete graph up to the substitutions N ^ N — 1 and t — 1 —>■ t needed due 
to different step rules and time origins. It also agrees with the conjecture formulated earlier in [25]. 
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in the scaling limit. is a generalized Gumbel distribution [38,39] which gives the 

standard type-I Gumbel distribution [33,34] 


T'l 


r ) = exp 


7) 


^-(t'+7) 


(4.42) 


when u = 0, i.e., at total covering (see hgure 4(c)). The discrete evolution of the 
probability distribution as total covering is approached {u = n — v = 0{1)) is shown in 
hgure 6. 


4-5. Crossover from Gumbel to Gauss 


In this section we study the crossover of the probability distribution (4.41) to Gaussian 
behaviour, when the deviation from total covering becomes larger, with n 3> m 3> 1. 
The generalized Gumbel distribution can be written as: 

= exp [/«+i(r' + 7 - Hu)] , fa{b) = -ab - e~’’ - Inr(a). (4.43) 

Let us expand fa{b) to second order in the vicinity of its maximum at bo: 

+ (4.44) 

One has: 


/:(6) = -a + e-^ 6o = -lna, ff{b) = -e-\ ff{bo) = -a. (4.45) 
When a ^ 1 the logarithm of the Gamma function has the following expansion [40]: 


In r(a) ~ a In a — a + - In 


27r\ 

a J 


(4.46) 


Finally: 

/a(&) ^-^In - |(6 + lna)^ (4.47) 

With 5 = r' + 7 — Hu, a = u + l^um (4.43), using (4.17) one obtains: 

5 + In a ~ r' + 7 — + In M ~ r'. (4.48) 


Thus, according to (4.47), the generalized Gumbel distribution in (4.43) can be 
approximated by 




g-UT'^/2 


(4.49) 


when n S> M ^ 1 in the vicinity of r' = 0. In order to compare to the Gaussian 
distribution %{T,f,x) = £3(t, y) in (4.32), obtained for x < 1 in the scaling limit, one 
has to change the scaling variable from t' in (4.37) to r in(4.28) so that: 


t — tisf^n{v) / _ 'T 

N “^ “ 7 ^' 


(4.50) 


When N is large, the behaviour of the probability distribution as a function of r is 
governed by the immediate vicinity of t' = 0. This justihes a posteriori the second 
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order expansion used above. The change of variable leads to the following Gaussian 
distribution: 


T 






u 


■ exp 


UT 


(4.51) 


y/N V 27riV 

This expression has to be compared to (4.32) when x is close to 1. In this limit, the 

logarithmic contribution to x{f) can be neglected so that 

I N N , , 

X(/, x) ~ —-- =-= — , 4.52 

f[l — X) n — V u 

and a perfect agreement between the two distributions is obtained. 


5. Discussion and outlook 

As mentioned in I, finite-size scaling results obtained for the fully-connected lattice 
(d = oo) are expected to be representative of the behaviour on periodic lattices above 
the critical dimension dc = 2. 

The effect of the restriction to a sublattice with n sites on the mean number of 
distinct sites visited sjsi^nif) given in (3.4) is quite simple. When normalized by the mean 
value of the total number of visits of the sublattice, tn/N, the ratio is independent of 
n and thus the same as for the full lattice. This property remains valid for periodic 
lattices in d = 1 to 3 as shown in [2]. 

A generalized covering time for random walks on graphs, the marking time, has 
been introduced in [42]. The walker marks a visited site i with probability Pi and the 
marking time is the time needed to mark all the sites. On the fully connected lattice, 
when the marking probability is uniform and equal to p, it is easy to verify that N 
has simply to be changed into N/p in the generating function (2.4). Thus, according 
to (3.15), the total marking time is given by NHn/p. 

The probability distribution for the number of distinct sites s visited up to time t 
in the subset (record number at t) and the probability distribution for the record value 
at a given record time lead to the same Gaussian density in the scaling limit. This is the 
probability density obtained in the mean-held approximation when correlations between 
the visits of different sites are neglected (see appendix B). But there is a remnant of the 
correlation effects in the variance which differs from the mean-held result by a term of 
order e”^”'. We have also obtained a Gaussian density for the record times at partial 
covering in the scaling limit. Finite-size corrections to the partial diherential equations 
leading to these results are of relative order Consistently, the asymmetry of 

the corresponding discrete distributions, as measured by the skewness, also vanishes as 
Ar-i/2, 

The record times are given by the sum of v independent, diherently distributed 
random variables, which are the lifetimes of previous records. Their variance scales 
as N at partial covering and as at almost total and total covering. Hence in the 
scaling variables, r for x < 1 and t' for x = 1, the time has to be divided by and 
N, respectively. The probability distribution of the record times is Gaussian at partial 
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covering, as expected since the central limit theorem applies for a sum of independent 
random variables with hnite variances. The generalized Gumbel distribution at almost 
total covering and the type-I Gumbel distribution at total covering are linked to the 
divergence of the variance of the lifetimes as s/n = x' —?■ 1 (see (4.2)). The main 
contribution to the record time then comes from lifetimes of records with number s 
close to n, with an exponential distribution given in (4.3). Note that this is not a 
standard case for Gumbel statistics which is usually associated with the distribution of 
extremes in a collection of random variables. 

The generalized Gumbel distribution, indexed by the deviation u = n — v from total 
covering, is crossing over to the Gaussian distribution obtained at partial covering when 
u increases. At hnite size, the crossover between the two regimes can be observed on 
the mean value, the variance and the skewness in hgure 5. The skewness in hgure 5(c) 
increases when v goes to n, reaching a value which converges from above to that obtained 
aX V = n m the scaling limit. 

The time evolution of s, the number of distinct sites visited by the walker on 
the sublattice with n sites, may be reinterpreted in different ways. Let us mention the 
example of a directed random walk with waiting times in ID where s in hgure 2(a) gives 
the position of the walker on the segment [0,n] as a function of time. The walker at s 
either takes a step forward with probability p{s) = {n — s)/N or waits with probability 
1 — p{s). The waiting times correspond to the lifetimes of the records. They are 
distributed according to JN,n{s, 1) in (2.6) and depend on the position of the walker. 
The mean velocity at s is p{s) so that the walker is slowing down as s increases and 
stops at s = n. In this interpretation the partial covering time is the hrst-passage time 
at an intermediate point s = v < n whereas the total covering time is the time of arrival 
at s = n. 

The random walk problem on the fully-connected lattice has a Russian dolls 
generalization with v walkers. The hrst walker performs a random walk on the fully 
connected lattice with N sites. The second is only allowed to take a step on the sites 
previously visited by the hrst, and so on and so forth. We are currently studying the ID 
transcription with a queue involving v random walkers in interaction on the segment 
[0,A^]. The walks are directed with waiting times depending of the distance between 
successive walkers. A walker is slowing down when approaching the preceding one in 
such a way that the initial order is always preserved §. 
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Note added in proof. It has recently been shown [45] that the generalized Gumbel 
distribution in (4.41) is universal for Markovian, non-compact random walks. 

Appendix A. Bivariate generating function for SN,n{s,t) 

Let us define a bivariate generating function SN,n{y,z) which is ordinary in y and 
exponential in 

= = . (A.1) 


t\ 1-^ — 

S=0 t=0 


s=0 t=0 

A Stirling number of the second kind can be expressed as (see I, equation (2.11)) 




(A,2) 


where A is the forward-difference operator such that A/(r 7 ) = /(?7 + 1) — f{ri). Making 
use of the definition (2.8) for a r-Stirling number of the second kind, one obtains 


r + t 
r + s 




k=0 


s! ' \k 

k=0 


rj^rt ^ = —A^(ri + rY 

s 

rt=0 


(A.3) 


T7=0 


Thus the generating function may be rewritten as: 

{z/Nf 


n c c oo 

n-y 




s=0 


_ ^z{l-n/N) 


E 

s=0 


t=0 

n 


t\ 


(t] + N — nf 


ri=0 




^s^zt]/N 


77=0 


(A.4) 


Since the result of the action of A on e^^/^ is a multiplication by (e^^^ — 1), one finally 
obtains: 

^ Q -!)]■= [l + - 1)]” , (A.5) 


When n = N equation (2.13) of I is recovered. 


Appendix B. Mean-field approximation for SN,n{s,t) 

The probability that a given site has never been visited up to time t is given by: 


(InY) = 


N-1 
N 


(B.l) 


Hence in the mean-held approximation, neglecting correlations between the visits of 
different sites, the probability distribution for the number of sites visited on the subset 
up to times t is the binomial distribution: 


SknM=(fj ll-toWl’toW"-* 


(B.2) 
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The mean value n[l — qN{t)], which is a sum of single-site terms, is in agreement with 
the exact result (3.4). The binomial distribution gives an extensive expression for the 


variance 


(B.3) 


= nqNit) [1 - qN{t)] , 

which, due to correlation effects, differs from the exact result in (3.6). 

In the scaling limit, the rescaled binomial distribution, ^{s,t), leads to a 

Gaussian density in the scaling variable a dehned in (4.10): 

,-aV(2<A') 


^(a,0') = 


\/27r0' 


(ffw) = e "'(1 


(B.4) 


The variance (j)'{w) is the rescaled variance As'%j^{t)/n in the scaling limit. It differs 
from (j){f,w) in (4.5) by a term proportional to /, of the second order in e~'^, which 
is a remnant of the correlations. This /-dependent correction to the mean-held result 
becomes less and less important as w increases, as shown in hgure 4(a). 


Appendix C. Partial differential equation for the probability density & 

We start from the master equation (2.12). Using the scaling variables in (4.1) and (4.10), 
the prefactors take the following forms: 


N -n + s s.i. , „ 

= 1 - /e 




n 


s -|- 1 s.l. p - 

= /e 


/VV 1 

iVi/2 + N ■ 


(C.l) 


N - . - . ^1/2 , ^ 

In the scaling limit n^ASj^ j^(^s,t) gives the probability density &[a{s,t), f,w{t)] which 
depends on s and t through a and w. A Taylor expansion in the variables s and t can 
be used on the left-hand side of the master equation (2.12) leading to: 

AV2 J dt 2 dt^ 

dG dG Id'^G Id^G d‘^G\ 
ds dt 2 ds'^ 2 dsdt) 


e = 1 - /e- 


/e- + 

J J^l/2 + JY 


6 - 


dt 2 ds"^ 

The needed partial derivatives are easily obtained and read: 
dG _ 1 dG dG _ 1 dG qq 92 q 

ds da ^ dt N dw 

d^G , o(iV-3/2) , ^^<3 


(C.2) 


1 d'^G 


da ’ 

-w q2q 


ds"^ Nf da"^ ’ 

+ 0(A-3/2^ . 


(C.3) 


dt"^ N da"^ ^ dsdt N da"^ 

Note that the second order expansion is sufficient to keep terms up to order N~^. 
Collecting coefficients of the different powers of (C.2), the leading non-vanishing 

contribution is of order N~^ and gives the partial differential equation (4.11). 


Appendix D. Partial differential equation for the probability density T 

With (4.1) and (4.28), the prefactors in the master equation (4.27) can be rewritten as: 
N — n + V 
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n-r + 1 s.l. . .1 


(D.l) 
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The scaling limit of t) is the probability density T[r(t, v), /, x{v)] depending 

on t and v through r and x. In the continuum limit, a Taylor expansion in t and v of 
the probability density on the left-hand side of (4.27) leads to: 


T = 


i-ni-n-jj 


f{l-x) + 
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^ 0% 1 d‘^%\ 

S —-1-I 

dt 2 dt^ J 

/ d% d% 1 dH 

T-^- 

\ dv dt 2 dv^ 


2 dt"^ dvdt) 


The partial derivatives have the following expressions: 
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1 
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(D.2) 


(D.3) 


Higher derivatives are of higher order in and can be ignored in (D.2). The 

leading non-vanishing contribution is of order N~^ and gives the partial differential 
equation (4.29). 


Appendix E. Expression of AfTjv j^(n — u,t) 


Inserting the explicit expression of the Stirling numbers of the second kind 
(equation (2.9) in I) 




into (2.8), one obtains: 
[r -|- mj ^ ml ^ 
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3=0 ' ^ k=0 

With t) in (2.10) one has the correspondence l = t — l,m = v — l = n — u — 1, 

r = N — n, so that: 

, n—u—1 


NTN,n{n -u,t) = 


n\ 
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^ ' j=0 


EM) 


n-u- 3-1 fn-U-l 


n-j 
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,(E.3) 


Finally, changing the sum over j into a sum over l = n — u — j — 1, equation (4.38) is 
obtained. 
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